Local superfluid densities probed via current-induced superconducting phase gradients 



David S. Hopkins, David Pekker, Tzu-Chieh Wei, Paul M. Goldbart, and Alexey Bezryadin 
Department of Physics and Frederick Seitz Materials Research Laboratory, 
University of Illinois at Urbana- Champaign, 1110 West Green Street, Urbana, Illinois 61801-3080, USA 

(Dated: July 10, 2007) 

We have developed a superconducting phase gradiometer consisting of two parallel DNA- 
templated nanowires connecting two thin-film leads. We have ramped the cross current flowing 
perpendicular to the nanowires, and observed oscillations in the lead-to-lead resistance due to cross- 
current-induced phase differences. By using this gradiometer we have measured the temperature 
and magnetic field dependence of the superfluid density and observed an amplification of phase 
gradients caused by elastic vortex displacements. We examine our data in light of Miller-Bardeen 
theory of dirty superconductors and a microscale version of Campbell's model of field penetration. 

PACS numbers: 74.50.+r, 74.78.-w 



Gradients in the phase of the superconducting order pa- 
rameter can be generated by external magnetic fields and 
sensed using a superconducting quantum interference de- 
vice (SQUID), invented by Jaklevic et al. [lj. Applica- 
tions of SQUIDS are numerous Q. As demonstrated by 
Clarke and others, SQUIDS may be used as magnetome- 
ters to detect remarkably small magnetic fields, e.g., mag- 
netic fields occuring in living organisms A voltmeter 
based on the SQUID, also developed by Clarke, is capa- 
ble of measuring 10 fV, i.e. it is 10 5 times more sensitive 
than a conventional semiconductor voltmeter Q . Micro- 
and nano-scale realizations of SQUIDs have been fabri- 
cated, e.g., in shunted Nb nanojunctions having a sub- 
micron hole [B| and, more recently, using nanotubes [f|. 
SQUIDs have also been used to study macroscopic quan- 
tum phenomena and devices, such as magnetization tun- 
neling 0, [H and phase qubits [§] . 

In this Letter we report on a molecular-templated ver- 
sion of the Jaklevic et al. phase-sensing experiment, in 
which DNA-templated superconducting nanowires are 
used instead of Josephson junctions to make a nanowire 
version of a SQUID (i.e. an N-SQUID). Our main accom- 
plishments are as follows. (1) By measuring the period of 
the resistance oscillations driven by an externally injected 
supercurrent, we have performed a direct, local measure- 
ment of the superfluid density and its dependence on 
temperature in the vortex-free regime. Our results can- 
not be consistently described by the phenomenological 
Gorter-Casimir and Ginzburg-Landau models. However, 
we do find good agreement with Miller-Bardeen's dirty 
limit of the microscopic BCS theory (loj , which was pre- 
viously examined for MoGe via a different method by the 
Lemberger group [lH . (2) We have observed and investi- 
gated the amplification of the current-induced phase gra- 
dient that occurs when pinned vortices are present in the 
cross-current-carrying lead. This amplification is brought 
about by a Lorentz force, which acts on randomly pinned 
vortices and results in their reversible displacement. In 
macroscopic settings, this physical phenomenon is the 
origin of the Campbell law, i.e. the dependence of the 
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FIG. 1: (color online) Schematic of the DNA-templated two- 
nanowire device. Two strands of DNA are stretched across 
a trench etched into SiN/Si0 2 on a Si chip. The molecules 
and banks are sputter-coated with Mo79Ge2i. A bias current 
Ib is applied through the nanowires and a cross current Jx is 
applied through the horseshoe lead, thus passing the contact 
points of the nanowires. A magnetic field B is applied per- 
pendicular to the plane of the leads and nanowires. One of 
the two metal-coated DNA molecules is shown on the SEM 
micrograph. 
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I15| . Our measurements provide a 
verification of the physics behind the Campbell law, but 
now at the microscale, as well as the capability of ob- 
taining the Labusch parameter, i.e. the average stiffness 
of the vortex pinning potential. 

Figure [T] shows the fabricated device. First, two MoGe 
nanowires are created by metal-coating a pair of sus- 
pended DNA molecules [16| . Then, a focused ion beam is 
used to cut one of the leads into a horseshoe shape, and 
thus define a narrow strip through which a cross current 
(Jx) can be applied. Independently, a bias current (Ib) 
can be applied through the nanowires, via which the lead- 
to-lead resistance is measured [TtJ • The typical width of 
the cross-current strip is of the order of 700 nm (e.g., 
633 nm for sample # 2), which is much smaller than the 
estimated zero-temperature perpendicular magnetic field 
penetration depth Aj_ = 46 /im. Thus, the analyses of the 
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FIG. 2: Resistance vs. temperature data for samples #2 (cir- 
cles) and #4 (triangles) . The open shapes correspond to zero 
cross-current measurements and the closed shapes correspond 
to a cross current that gives a resistance close to a maximum 
(3.8 for sample #2 and 5.9 /xA for sample #4). The fits 
are obtained with the following fit parameters: sample #2 - 
Rn = 800 Q, Jd(0) = 376 nA, J c2 (0) = 151 nA, T cl = 2.45 K, 
T c2 = 1.86 K; and sample #4 - R N = 700 fi, I c i (0) = 1066 nA, 
/ c2 (0) = 317nA, Td = 3.04K, T c2 = 2.55K jSJ. 



results will be carried out under the assumption that the 
cross current flows uniformly in the cross-current strip. 
The width of the larger of the two leads, 21 = 17.33 fim 
is smaller than Aj_, and thus any external field will pen- 
etrate both leads without significant attenuation. More- 
over, we estimate that the effects of the magnetic field 
generated by the cross current are negligible. The length 
of the nanowires, denoted by 6, is 100 nm for sample # 
2, and the distance between them, denoted by 2a, which 
determines the spatial resolution of the device, is 3.13 /im 
for this sample. Devices with a better spatial resolution 
(e.g. with 2a 265 nm) were demonstrated in a different 
setup 0. 

At 7x = we observe a broad resistive transition 
in the nanowires, as the temperature T is decreased. 
Moreover, as Ix is increased, this transition periodically 
broadens and narrows back to its Ix = breadth. We 
show this transition in Fig. [2] for two samples, each mea- 
sured at two distinct values of Ix that correspond al- 
most exactly to the minimum and maximum observed 
transition breadths. The T dependence of the resis- 
tance follows the theory of thermally activated phase 
slips (TAPS) US [2(1, extended to the two-wire 
case [ill [2lj]. This extended theory involves a modi- 
fied free-energy barrier for phase slips, which accounts for 
the inter- wire coupling mediated through the leads. In 
the short- wire limit, the two- wire device has the current- 
phase relation of a single Josephson junction with an ef- 
fective critical current 
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FIG. 3: Resistance vs. cross-current data for sample #2 at 
temperatures ranging from l.OKto 1.8Kin0.1K increments. 
The solid lines are fits using the same nanowire parameters 
as the fit for the R-T data in Fig. [2] and for the period vs. T 
data in Fig. [4] The dashed vertical line corresponds to the 
position of the peak at 1 K, and is an aid to see the decrease 
in period at higher temperatures. The arrows point to peaks 
and show the change in the oscillation period with T. 



nanowires S (see Fig. []}. For the cross-current exper- 
iment, S = J" a V0 • dr, where 4> is the phase of the 
superconducting order parameter and the integral runs 
between the nanowire ends via the horseshoe lead. We 
note that the oscillation period, which is the main focus 
of this Letter, is independent of whether the nanowires 
are operating in the short wire limit or not. Moreover, 
J c i and J C 2 are the critical currents for the nanowires, and 
are given by J cli2 = 7 c i, 2 (0)[l - T/T cl , 2 ] 3 / 2 f^. We ob- 
tain the device resistance R from the damped Josephson 
junction formula 0, [24| (Eq. (9) in Ref. [3|), where we 
have used the normal-state resistance i?N of the two par- 
allel nanowires for the effective shunt resistance of the 
junction. Choosing <5 = and 5 = ir, we fit the lower 
and upper curves in Fig. [51 As shown in Fig. [31 the re- 
sistance oscillates as a function of Ix, having a period 
on the order of 7/xA and an amplitude that is maximal 
in the middle of the resistive transition. As the arrows 
in Fig. [3] indicate, the period of the resistance oscilla- 
tion is temperature-dependent, in contrast with the case 
of magnetic- field-induced oscillations, which appear to 
be temperature independent [lij ]. Thus, we have direct 
proof that the oscillation is not due to a magnetic field 
induced by the cross current. 

The period of the resistance oscillation is determined 
by the condition that the phase gain S be an integer 
multiple of 27T. For the case of a uniform cross cur- 
rent in a thin-film strip of width w and thickness d, 
is related to Ix through the superfluid density n s : 
= 6 /2a = (Ix/n s )(2m/wdeh) 1 where m and e are 
the electron mass and charge. Thus, by measuring the 
period of the resistance oscillation vs. T, we obtain the 
superfluid density in the strip carrying the cross current 



3 



0.2 



0.4 



0.6 



.9 

0.8 
1 

0.8 
.6 




fat 



» Sam pit #2 

\ xSample#4 
\+ Sample #6 




0.2 0.4 
T 



0.6 



FIG. 4: (color online) Normalized superfluid density n a 
vs. T/T c . The solid line represents Miller-Bardeen's theory 
(Eq. [2}; the dotted line is Gorter-Casimir (GC) model; the 
dotted-dashed line is Ginzburg-Landau (GL) formula. Top: 
The n s (0) was chosen to optimize Miller-Bardeen's fit, the 
respective scaling parameters A7x(0) being 7.39 n A (sample 
#2), 11.70 fiA (sample #4), and 7.30 /iA (sample #6). The 
values of T c used here, 5.64 K, 5.72 K, and 5.77 K for samples 
#2, #4, and #6, respectively, have been measured indepen- 
dently. Bottom: The same set of data is shown, but n s (0) 
was chosen to optimize the GC fit (top group of the points) 
or the GL fit (bottom group of points). 
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FIG. 5: (color online) Resistance vs. cross current measured 
at various values of the field at 1 K. For clarity, all the curves, 
except that of mT, have been offset vertically by an integer 
multiple of 170 Q. The 100 mT and 200 mT curves have 
been horizontally magnified by a factor of 4, since the pe- 
riod became rather small at these fields. The curves show 
multiple random jumps at higer fields. By measuring the re- 
sistance over a small range of cross currents several times, it 
was possible to measure at least half a period of oscillation, 
as indicated by the horizontal arrows. 
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The normalized period of the resistance oscillation (and 
hence the normalized superfluid density) vs. T is shown 
in Fig. 2J To make the fits we have used Miller-Bardeen's 
result, applicable to the dirty superconductor case [Io| . 
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where T c is the critical temperature, n s (0) is the zero- 
temperature superfluid density, and A(T/T C ) is the uni- 
versal BCS gap relation 25]. The fit to Miller-Bardeen's 
formula was performed by allowing n s (0) (and hence 
A7x(0)) to be an adjustable scaling parameter. It gives 
better agreement than do the alternative theories, such 
as Gorter-Casimir and Ginzburg-Landau. In fact, the fit 
to Miller-Bardeen's theory can be further improved by 
slightly varying T c . Almost perfect Miller-Bardeen-type 
fits are obtained with T c = 5.22, 5.94, and 5.56 K for sam- 
ples #2, #4, and #6, respectively. Combining the ob- 
tained fitting parameters Aix(0), I c i, T c \, I C 2, and T C 2, 
we produce the theoretical curves in Fig. without any 
additional adjustable parameters. 

Having shown that our device is capable of sensing re- 
sistance oscillations due to phase gradients created by 
external supercurrents, we now turn our attention to the 
application of magnetic fields > H c \, which induce vor- 
tices in the thin-film leads and thus create additional 



phase gradients. The penetration of vortices into meso- 
scopic samples was studied prev iously using other meth- 
ods (see, e.g. Refs. 27, H, 29|). All the measurements 
presented in the following were performed on sample 
#2. Figure [5] shows the resistance as a function of cross 
current at several values of applied magnetic field. At 
low fields, when vortices are absent, the cross-current- 
induced oscillations are phase shifted with respect to 
one another for distinct magnetic fields but with no de- 
tectable difference in the oscillation period (see curves 
labeled mT and 3.4 mT). This shift is due to the addi- 
tional phase gradient caused by the Meissner currents as- 
sociated with the applied magnetic field. At fields slightly 
above H c \ for the wider lead (i.e. not the cross-current 
lead), it is possible to see one of two types of resistance 
traces. Sometimes, we observe noisy oscillations with a 
smaller than expected amplitude (e.g. 4.5 mT) but with 
the same period as that at B = 0. We believe that this 
is due to one or more vortices that are rapidly wander- 
ing near one or both of the nanowires, i.e. a version of 
motional narrowing. At other times, we observe the type 
of oscillation as observed below H cl , but with occasional 
jumps (6.7 mT) which become more frequent when the 
field is larger. At fields above H c i for the cross-current 
carrying lead (i.e. the lead with the larger H c i), the pe- 
riod begins to decrease, as shown by the highlighted seg- 
ments for the 100 mT and 200 mT traces in Fig. [5] As 
vortex jumps are prevalent at these fields, and our goal 
is to measure the period for a given vortex configuration, 
we determine the period by measuring it within "quiet 
regions" (i.e. current intervals exhibiting no jumps). Ex- 
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FIG. 6: (color online) Inverse of the oscillation period vs. field. 
The black line is from Campbell's formula with Ax = 46 fim 
and k — 370 N/m 2 . The inset shows the low-field data. 



amples of such regions are indicated in Fig. [5] by horizon- 
tal arrows for the 100 and 200 mT traces. 

We have investigated the resistance oscillation period 
as a function of magnetic field at four temperatures by 
making several scans of Ix in "quiet regions," in which it 
was possible to observe at least one peak and one val- 
ley several times before a vortex jump occurred. At 
temperatures below ~ 700 mK, at which the device has 
an undetectably small zero-bias resistance, we obtained 
the period from large bias-current resistance measure- 
ments (i.e. beyond the linear response regime). At each 
field and temperature, several measurements of the pe- 
riod were recorded and averaged. Smoothed curves of 
the inverse period, obtained via a moving average, are 
shown as functions of magnetic field in Fig. [5] We ob- 
serve a roughly linear increase of the inverse period; com- 
pared to its zero-field value, the period at B — 3 T is 
some thirty times smaller. Such behavior is analogous 
to the linear dependence of the perpendicular penetra- 
tion depth Aj_ on applied magnetic field (in the regime 
in which the vortex density is pro portional to the field), 
i.e. the Campbell formula [l2l Il3l. adapted for thin 
films: A c ,_l = (A 2 + ($ £/MCn))/d. Here, k(T) is the 
vortex pinning force constant (i.e. the Labusch parame- 
ter [HI). In the present setting, the phase fields from the 
vortices, displaced due to Ix, enhance the phase gradi- 
ents resulting from the cross current, resulting in a period 
for resistance oscillations that shortens with B. We can 
therefore tie the Campbell prediction to the measured pe- 
riod by assuming the usual relation n s — m/ ^oe 2 d\c,±- 
From this and Eq. ([1]) , we expect that the inverse period 
should increase linearly with magnetic field for our thin- 
film leads, (ATx) -1 = (2a/wd) (a* A 2 /$o + (B/k{T))). 
We estimate that k « 370 N/m 2 ; other materials have 
yielded, e.g., > 100 N/m 2 and - 10 5 N/m 2 [30]. Fur- 
thermore, the inferred value of X±(B = 0) is consistent 
with an estimate made using bulk penetration depth and 
film thickness data [3l[. As Fig. [5] shows, the inverse pe- 



riod does indeed increase linearly, except at higher tem- 
peratures and fields. In those regimes, inter-vortex in- 
teractions become important and hence, k can be effec- 
tively reduced. At high temperatures, vortices explore a 
larger area and are less sharply pinned. For fields below 
H c i ~ 0.02 T, no vortices exist in the horseshoe lead, and 
the period does not depend on the field. The deviation 
from linear ficld-dcpcndcncc observed slightly above H c \ 
(Fig. [6] inset) arises because in this regime the vortex 
density grows sublinearly with the field. 

Regarding the future applications of our N-SQUID gra- 
diometer, we remark that it can be used not only to mea- 
sure the local superfluid density, e.g., in connection with 
the fate of the Josephson effect in ultra-thin film sys- 
tems [HI, but also to shed light on settings in which lo- 
cal fluctuations of the superfluid density have been pre- 
dicted, such as highly disordered s-wave 2D supercon- 
ducto rs | 33I| and superconductors having magnetic impu- 
rities [34j. 
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